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ABSTRACT 

HH SON is one of the Herbig-Haro objects that have associated quiescent dense clumps. We report 
CO and CS BIMA observations that reveal star formation within the HH SON dense clump. The CO 
emission reveals clearly a bipolar molecular outflow centered on the dense clump. The CS emission traces 
a ring-like structure of radius ~ 0.24 pc. The CS kinematics shows that the ring is collapsing with an 
infall speed of ~ 0.6 km s~^. The required mass to produce the collapse is in agreement with previous 
ammonia observations of the 20 Mq core, which is embedded within the CS structure. However, we 
cannot discard that the ring structure is expanding driven by protostellar winds, if the CS abundance 
if unusually high and the CO momentum rate is much higher than that measured, due to inclination 
and optical depth effects. The properties of the molecular outflow and of the dense core suggest that it 
harbors a Class object. There are also signatures of interaction of the HH SO/Sl/SON outflow with the 
dense gas. In particular it is possible that the HH SO/Sl/SON outflow has triggered or at least speed up 
the star formation in this region. 

Subject headings: ISM: individual (HH SON) — ISM: jets and outflows — ISM: kinematics and 
dynamics — stars: formation 



1. INTRODUCTION 

HH SO/Sl/SON is a spectacular and very luminous HH 
outflow spanning about 5 pc and exhibiting proper mo- 
tions, close to the exciting source, of up to 1400 km s^^ 
(Reipurth & Graham 19SS; Martf, Rodriguez, & Reipurth 
1993, 1995, 1998; Heathcote, Reipurth & Raga 199S; Moli- 
nari, Noriega- Crespo & SpinogHo 2001). HH SON, the 
northern counterpart of HH SO/SI, is optically invisible, 
due to the high extinction at this position (Martf et al. 
1993), but its radio and far- infrared properties confirm its 
HH nature (Martf et al. 1993; MoHnari et al. 2001). Far- 
infrared ISO observations show that the strong FUV field 
radiated by the HH 80/81/80N outflow have induced the 
formation of a PDR (photo-dissociation region) in its im- 
mediately surrounding region (Molinari et al. 2001). 

Molecular line observations towards HH SON (Girart et 
al. 1994; Girart, Ho & EstaleUa 199S) showed that it be- 
longs to the class of Herbig-Haro (HH hereafter) objects, 
typically of high excitation, that are found to have starless 
quiescent dense clumps ahead of them. The strong radi- 
ation arising from the HH objects penetrates the clumps, 
releasing the icy mantles from the dust grains and induc- 
ing a photon-dominated chemistry, which causes signifi- 
cantly altered chemistry (Taylor & Williams 1996; Viti & 
Williams 1999). Recent observations towards HH 2 have 
confirmed this chemical scenario (Girart et al. in prepara- 
tion). However, the presence of IRAS 1S163-2042 associ- 
ated with the HH SON dense clump suggested that it could 
be tracing an embedded star (Girart et al. 199S). 

This paper presents CS and CO BIMA observations that 
are part of a molecular line survey towards HH SON whose 



complete results will be reported in a forthcoming paper. 
The observations presented in this paper show infall and 
outflow signatures associated with the dense clump down- 
stream of HH SON. The high angular resolution provided 
by millimeter interferometers, as in this paper, have al- 
lowed astronomers to resolve spatially infall kinematical 
signatures in some low mass star forming cores (Zhou, 
Evans & Wang 1996; Ohashi et al. 1997; Saito et al. 1996; 
Momose et al. 1998; Hogerheijde 2001). Typically, the in- 
falling gas arises from scales of ~ 1000-3000 AU and the 
mass infall rates are 1-10x10"^ Mq yr~^. Interestingly, 
the infall signatures found in HH SON are similar to that 
from the dense cores in L1544 (Ohashi et al. 1999) and 
HH 1/2 (Torrelles et al. 1994), which show infall signa- 
tures in the form of an almost ring-like structure in the 
position- velocity diagram on even larger scales, 1.5 x 10* 
and 4.5 x 10* AU, respectively. 



2. OBSERVATIONS AND RESULTS 

The observations were carried out between 1999 May 
and September with the 10-antenna BIMA array^ at the 
Hat Creek Radio Observatory. The phase calibrator used 
was QSO 1733-130. Absolute flux calibration was per- 
formed by observing Uranus and Mars. The phase cen- 
ter of the observations was a(J2000) = 1S'^19'"18!62; 
(S(J2000) = -20°40'55". Two different frequency setups 
were used in order to observe the CS J=2-l and CO 
J=l-0 lines, which provided a spectral resolution of 0.3 
and 0.5 km s^*, respectively. Maps were made with the 
(u, v) data weighted by the associated system tempera- 
tures. 



^ Departament d'Astronomia i Meteorologia, Universitat de Barcelona, Av. Diagonal 647, 08028 Barcelona, Catalunya, Spain; jgirart@am.ub.es 

2 Department of Astronomy, University of Illinois, 1002 W. Green St, Urbana, IL 61801, USA 

^ Department of Physics and Astronomy, University College London, London, WCIE 6BT, England 

* Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

^ The BIMA array is operated by the Berkeley-IUinois-Maryland Association with support from the National Science Foundation. 



2 



Star Formation Signatures in the Condensation Downstream of HH SON 



2.1. The CS J =2-1 Emission 

The CS emission arises from an elongated structure in 
the NW to SE direction {PA ~ 125°), consisting of two 
well differentiated structures morphologically and kine- 
matically (see Figures 1,2): (1) The main core is centered 
at the position of the ammonia core (Girart et al. 1994). 
The position-velocity plot along the major axis of the CS 
emission shows an ring-like structure, roughly symmetri- 
cal about Vlsr — 11-4 km s~^. This structure can be 
explained by inward (or outwards) motions in the conden- 
sation (sec § 3.1) and resembles that from the starless core 
L1544 (Ohashi et al. 1999). The physical scale of the dense 
main core traced by the CS emission, ~ 0.48 x 0.17 pc (as- 
suming a distance of 1.7 kpc: Marti et al. 1993) is larger 
than that from the NH3 (1,1) (Girart et al. 1994) and the 
HCO+ J=3-2 (Girart et al. 1998). This could be due to 
the higher critical density of these two lines with respect 
to the CS J=2-l transition. (2) A small redshifted "arm" 
(hereafter RSE arm) is located SE of the main core and 
~ 15" north of the 6 cm continuum emission. The RSE 
arm c;xhibits a clear linear velocity gradient with increasing 
velocity towards the east, i. e., as it gets closer to HH SON. 
Table 1 shows the physical parameters of the core traced 
by the CS J=2-l line. Of the total mass measured, 12 Mq 
(assuming X(CS) = 2 x 10~^), around 10 M© belong to 
the main core, and 2 Mq to the RSE arm. 

2.2. The molecular outflow 

The blueshifted and redshifted CO emission traces a 
clearly bipolar outflow along the E-W direction (lobes BE 
and RW) and approximately centered at the peak inten- 
sity of the CS integrated emission (Fig. 1). The presence 
of molecular outflows is a clear signpost of star formation 
within the dense core. The outflow axis is not perpendic- 
ular to the major axis of the main core. However, our an- 
gular resolution traces scales (^ 2 x 10^ AU) much larger 
than the scales at which the coUimation process occurs, 
~ 1 AU (Shu et al. 1995). Table 1 shows the physical pa- 
rameters of the molecular outflow. There is marginal evi- 
dence of another bipolar outflow in the N-S direction and 
also centered in the dense core. In addition to the molecu- 
lar outflow, there is a redshifted clump that coincides not 
only spatially with the RSE arm but also kinematically, 
suggesting that this CO emission traces the same gas as 
the CS at this location. 

3. DISCUSSION 

3.1. The collapse of the main core 

The position-velocity plot along the major and the mi- 
nor axes of the CS emission (Fig. 2) clearly suggests that 
the gas is either collapsing or expanding radially and that 
the emission arises from a ring-like or toroidal structure. A 
similar structure has being observed towards the starless 
core L1554 from CCS observations (Ohashi et al. 1999). 
The ratio between the major and minor axes of the CS 
emission gives an inclination angle of > 70°. As pointed 
out by Ohashi et al. (1999), the fact that the emission 
is detected towards a ring-like structure is probably due 
to a chemical effect, i. e., a decrease of the CS abundance 
toward the center, possibly because of gas depletion onto 
grains at higher densities (Rawlings et al. 1992). 



We considered a model of a spatially thin ring with both 
infall and rotation, similar to that of Ohashi et al. (1997, 
1999), seen edge-on by the observer. The parameters of 
the model were the inner and outer radii of the ring, i?i„n 
and i?out, the infall and rotation velocities, l^niaii and 
\4ot) and the intrinsic line width AVcore- The position- 
velocity diagram was computed along the projected major 
and minor axes of the ring, assuming that the intensity 
was constant for i?inn < R < Rout, to compare with the 
observed kinematics. The best fit was obtained for veloc- 
ities Vinfaii = 0.59 km s~^, K-ot ^ 0.20 km s~^, and radii 
of the ring i?inn = 0.19 pc, Rout = 0.29 pc (see Table 1). 
The observed and synthetic position- velocity diagrams are 
shown in Fig. 2. This model also applies if the ring is ex- 
panding instead of contracting. In this case, the (ierivc;d 
infall velocity should be taken as the expantion velocity. 

Can the protostellar wind cause the expansion of a ring- 
like structure? The momentum rate of the ring (see Ta- 
ble 1) is 2 X 10"^ Mq km s"^ yr"^ On the other 
hand, the momentum rate of the protostellar wind, assum- 
ing that it is the same as that from the molecular outflow, 
is ~ 1 X 10~^ Mq km s~^ yr~^. Although both values are 
similar, only a small fraction of the wind momentum rate 
would be transfered to the dense core. First, the dense 
core, as traced by the CS, shows a relatively flattened 
structure. Assuming that the ring structure traced by the 
CS is a torus, then the inner and outer radius derived from 
the model will imply a semi-aperture angle of the torus as 
seen from the center of 12°. Second, the momentum rate 
distribution of the wind is expected to be anisotropic, i. e., 
most of the momentum rate of the wind is carried out 
along the axis of the outflow, e.g., Shu et al. (2000). The 
X-wind model (Shu et al. 1995) predicts that at a distance 
of 0.24 pc the fraction of wind mass moving at equatorial 
angles lower than 12° is about 2% of the total wind mass. 
We expect the momentum rate to scale similarly. Thus the 
wind that could interact with the dense core has roughly a 
momentum rate 2 x 10"'' Mq km s~^ yr~^, which is much 
lower than that measured from the CS kinematics, assum- 
ing a typical CS abundance of X(CS) = 2 x 10~^. Thus, 
the wind would be able to cause the expansion of the ring 
only if the CS has an unusually high abundance, and if the 
true momentum rate is much higher than that measured, 
because of the outflow inclination and the optical depth. 

Can the ring-like structure be in contraction? Girart et 
al. (1994) found that the ammonia emission arises from 
a smaller region than the CS J=2-l and traces a 20 Mq 
core (see § 2). This mass is similar to that required to 
cause free-fall collapse for the infall velocity measured at a 
radius of ^ 0.24 pc. However, the gas pressure and mag- 
neticc field should be taken into account. For the intrinsic 
line width obtained from the model, the virial mass with- 
out magnetic field for the CS main core is 9 Mq. This 
indicates that in the absence of the magnetic field the gas 
pressure cannot prevent the collapse. Assuming that the 
ring-like structure is contracting, we can provide a rough 
upper limit on |B|, i.e., the main core should be super- 
critical and the total virial mass (including the magnetic 
fields) should be no more than the free-fall mass. Follow- 
ing the formulae given by Crutcher (1999) and McKee et 
al. (1993) we find that for a mass of 20 Mq the collapse is 
possible for |B| < 50/xG. Magnetic fields of < 50/xG have 
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been measured in several molecular clouds with similar 
sizes and line widths as those of HH SON (Crutcher 1999). 
Thus, it seems that there is enough mass to cause collapse 
while, on the other hand, an expansion would require an 
unusually high CS abundacc and a high correction to the 
observed CO momentum rate. Therefore, we favour that 
the CS kinematics is tracing infall. 

Compared with most of the low mass regions where a 
kinematical infall signature has been found (see § 1), the 
HH SON main core; has a larger mass infall rate (sec; Ta- 
ble 1). This is probably because it is also a more mas- 
sive region. Another example is the infalling NH3 core in 
HH 1/2 (Torrelles et al. 1994), where the infall material 
also comes from a ring-like structure (toroid) with a simi- 
lar size and mass as the HH SON main core. Interestingly, 
the HH SON and HH 1/2 star forming cores show clear 
differences with the other collapsing ring-like structures 
detected so far. Here, the infalling motions arc clearly 
supersonic with no rotation observed, whereas the infall 
motions in the starless core L1544 are subsonic and of the 
order of the rotation speeds (Tafalla et al. 199S; Ohashi 
et al. 1999; Williams et al. 1999). These differences could 
be due to an evolutionary effect (L1544 is still a starless 
core) and/or to the different masses involved in the col- 
lapse. The powering source of HH 1/2 is a protostar with 
an far-IR luminosity of ~ 50 (Harvey et al. 1986), but 
what about HH SON? 

3.2. The evolutionary stage of the HH SON main core 

IRAS 18163-2042 has a luminosity of ~ 20 Lq (Girart 
et al. 1994), but because of the large position uncertainty 
it is possible that both the PDR region heated by HH SON 
(Molinari et al. 2001) and the protostar embedded in the 
infalling core contribute to this Imninosity. Therefore, this 
value is a rough upper limit of the protostellar luminosity. 

A detailed study of the outflow activity, from a homoge- 
neous sample of low-luminosity embc^dded YSOs, done by 
Bontemps et al. (1996) found that the relation between the 
outflow momentum flux, the envelope mass and the bolo- 
metric luminosity depends on the evolutionary stage of the 
protostars. In particular, the relation between Fco c/^boi 
and Mcnv/^^bol' which should be almost free of any dis- 
tance and luminosity effects, is useful to distinguish the 
evolutionary stage of the low mass YSOs, that is, whether 
they are in the Class or I stage. The values of the mea- 
sured envelope mass and the corrected momentum rate, 
Fco = 10-P (the same inclination and opacity correction 
factor used by Bontemps et al. 1996), suggest that, for the 
possible range of bolometric luminosities (<20 Lq), the 
HH SON dense core is in the Class evolutionary stage. 
Bontemps et al. (1996) also found that the momentum flux 
and the envelope mass are well correlated for both Class 
and I objects. However, taking into account the momen- 
tum flux of the HH SON molecular outflow, the dense core 
is an order of magnitude more massive than what is ex- 
pected from this correlation. The excess of envelope mass 
with respect to the outflow strength could be explained if 
there is multiple star formation going on (e. g., if at unre- 



solved scales the core splits into several clumps, with some 
of them still pre-stellar) . 

3.3. Interaction between HH SON and the molecular cloud 

There are two observational features that suggests that 
HH80/81/80N outflow may be interacting with the am- 
bient gas surrounding the infalling main core: First, the 
clear velocity gradient of the RSE arm seen in both CS 
and CO emission and its location just above HH SON (see 
Fig. 2). Second, the BE lobe of the molecular outflow 
spreads out to the north at the position where it coincides 
spatially with the axis of the HH80/81/80N outflow (see 
Fig. 1: the arrow line indicates the outflow direction of 
HH SON). This scenario, where HH SON and the dense gas 
are almost co-spatial, suggests that the infalling dense core 
and the RSE arm may still be affected chemically by the 
strong UV radiation of HH SON, as suggested by Girart 
et al. (1994, 1998), since this radiation is so strong that it 
has created a PDR region in its closest vicinity (Molinari 
et al. 2001). 

It is unlikely that HH SON triggered the star formation 
within the main core: assuming a HH SON propagation 
speed of 100-300 km s~^, and using the typical age of a 
Class source for the YSO embedded in the main core, 
~ 10'' yrs, HH SON would have to be about 1 to 3 pc 
away from the main core when the collapse was initiated. 
Yet, the material from the HHSO/Sl/SON flow was Hkely 
present near to the CS core before HH SON arrived to its 
proximity: a giant bow shock is observed dowstrcam of 
HHSO/Sl (Heathcote et al. 1998) and the analysis of the 
line emission from HHSO/Sl/SON shows that the shocks 
arise at the interface between two fast-moving flows (Moli- 
nari et al. 2001). Therfore, it is still possible that the 
HHSO/Sl/SON outflow has triggered or at least speed up 
the star formation. There are many uncertainties in our 
understanding of the dynamics of collapse; these dynamics 
are determined by the balance of gravity against pressure. 
Once the collapse is triggered, it is yet not clear how in- 
fall proceeds siiic;c; thermal and non-thermal pressure may 
support the clump from collapsing, or at least retard an, 
often assumed, free-fall collapse (Rawlings et al. 1992). In 
the case of HH SON, the UV radiation from the shock front 
may play an important role in the dynamics of the collapse 
of the dense core since it will highly affect its cooling func- 
tion. Viti & Williams (1999) showed that the radiation 
generated in the HH shock drives a rapid chemistry in the 
dense core, by lifting the icy mantles and therefore en- 
hancing the gas with hydrogenated species, such as NH3 
and H2O; the latter, for example, will undoubtedly cool 
the gas. This would 'encourage' the infall. This scenario 
will be investigated in a further paper. 
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Fig. 1. — Superposition of the: gray scale map of the VLA 6 cm continuum emission from Marti et al. (1993); grey thick contour map of the 
CS J=2-l integrated emission over the 9.8-13.6 km s~^ vi^sR interval; dotted thin contour map of the averaged CO J=l-0 blueshifted emission 
over the 4.7-9.7 km i'lsr interval; dashed thin contour map of the averaged CO J=l-0 redshifted emission over the 12.7-17.7 km 
^LSR interval. CS contours are 0.5, 1.0, 1.7, 2.4, 3.1 Jy beam~^ km s~^. CO contours axe 3, 5, 8, 7, 11 x 0.14 Jy beam~^, the noise level 
of the maps. The cross indicates the position of the ammonia peaJf intensity (Girart et al. 1994). The solid arrow line shows the direction 
of the HH 80/81/80N flow. The half-power contour of the CS and CO synthesized beams (15'.'6 x 7'.'1, PA = 3° and 14'.'9 x S'.'O, PA = 3°, 
respectively) are shown on the bottom left-hand corner, with the CS beam to the right. The top panels show the CO J=l-0 spectra averaged 
over the BE and RW lobes. 
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Fig. 2. — Bottom left panel: Position- Velocity (PV) plot of CS J=2— 1 emission along the major axis of the dense core, PA = 125°, with the 
position offset relative to «(J2000) = 18'' 19™ 18! 13; 5(J2000) = -20°40'40'.'4. Bottom right panel: PV plot of CS J=2-l emission along the 
minor axis of the dense core. Top panels: Synthetic emission of the PV plots along the major (right) and minor (left) axis of an infalling ring 
with an infall velocity of 0.59 km s"'^, no rotation, an inner and outer radius of 0.19 and 0.29 pc, respectively, and line width of 0.8 km s~^. 
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Table 1 

Physical Parameters of the HH SON region* 



12[X/2 X 10-9]-i 
11.43 ^ 
0.7^ 
0.19 ^ 
0.29 ^ 
< 0.20 b 
0.59 ^ 

2 X 10"^ 

3 X 10"^ 



0.11 
~0.4 
10 

3.6 X 10^ 
0.35 

1.0 X 10-5 

LI X K)--' 



Assuming optically thin emission, LTE and no cor- 
rection for inclination of the outflow 

''Prom the infalling ring model (see § 3.1). 



Molecular Dense Core: 
M{Mq) 

Vcore (km S ^) 
AVcore (kni S"^) 
-dinner (pc) 
^outor (pc) 
Vrotation (km S"^) 
Mnfall (km S"^) 

Anfaii (M© km s-1 yr-i) 
Mnfaii (Mq yr^i) 

Molecular Outflow: 

M (Mq) 
^max (pc) 
Vmax (km S^l) 

Td (yr) 

P {Mq km s-1) 

P (Mq km yr-i) 

Lnirrh (L . ) 



